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, ^ , Abstract 
, 

' We study the similarity renormalization scheme for hamiltonians to the fourth order in 

perturbation theory using a model hamiltonian for fermions coupled to bosons. We demon- 
strate that the free finite parts of counterterms can be chosen in such a way that the T-matrix 
is covariant up to the fourth order and the eigenvalue equation for the physical fermion re- 
duces to the Dirac equation. Through this choice, the systematic renormalization scheme 
I reproduces the model solution originally proposed by Glazek and Perry. 

o 

. 1 Introduction 

Our study of the similarity renormalization scheme for hamiltonians ^ is done in the model 
which consists of only two sectors in the Fock space. This great simplification of the space of 
states allows complete analysis of the renormalization scheme and still includes typical factors 
• and divergences that appear in quantum field theory. Our model is based on Yukawa theory. 

The hamiltonian of Yukawa theory truncated to one fermion and one fermion plus one 
k> ■ boson Fock sectors leads to infinities in the fermion-boson T-matrix. Therefore, we introduce a 

j_j ' cutoff A for the momentum transfer in the interaction part of the hamiltonian. The similarity 

transformation allows us to construct counterterms in the initial hamiltonian in such a way that 
the renormalized hamiltonian gives finite and cutoff independent results for the T-matrix. We 
construct renormalized hamiltonians using expansion in powers of the effective fermion-boson 
coupling constant and including terms up to the fourth order. 

In the similarity renormalization scheme, one constructs effective hamiltonians Hx which are 
functions of the width A. H\ is obtained from the initial hamiltonian H\ with the imposed cutoff 
A and added counterterms by a unitary transformation. The transformation and counterterms 
are found order by order in perturbation theory using the requirement that matrix elements of 
H\ are independent of the cutoff A when the cutoff goes to infinity. 

To find the unknown finite parts of the counterterms we calculate the T-matrix for fermion 
- boson scattering. The condition that the T-matrix is covariant can be satisfied and it implies 
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relations between the finite parts of different counterterms. We also demand, that the physical 
fermion is described by the Dirac equation with the fermion mass equal to the fermion mass 
term in the fermion-boson sector. This demand also provides a relation between the finite parts 
of counterterms and it is called the threshold condition [^. 

The model hamiltonian we study was originally considered by Glazek and Perry Q. They 
guessed the form of counterterms which remove divergences in T-matrix and they obtained 
covariant results for the T-matrix to all orders. 

Our main question about the model was if the systematic similarity calculation, carried out 
in perturbation theory, produces the same solution to the hamiltonian renormalization problem 
as guessed by Glazek and Perry. The cutoff in the model is limited by the triviality bound but 
one can assume that the coupling constant is small enough for reliable use of the perturbation 
theory. 

Section 2 presents the model. Sections 3 and 4 describe its renormalization to the fourth 
order. Section 5 explains connection with Ref. We conclude in Section 6 and Appendix 
contains key details. 

2 Model 

The initial hamiltonian is a light-front hamiltonian for Yukawa theory projected on two Fock- 
space sectors: one with a fermion and one with a fermion and a boson. 

HA = Hof + Hofb + HY + H+ + XA, (1) 

where the free part is 

[p,k]\pa,k){pa,k\i — — ^ k+ J ' 
boson creation and annihilation vertices are 

Hy = gYl f[pi,P2,mA^-M^^^f,)2i27rf6'{pi-p2-k)x 

0"l,0-2 

x[\p2a2,k){piai\u{p2,a2)u{pi,ai) + h.c] = + , (4) 

and the seagull term is 

H+ = g^Y. l[pi,P2,h,k2]e{A^-Ml)e{A^-Ml)2{27Tf6Hp2 + k2-pi-ki)x 



0"l,0'2 
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x\p2a2,q2)(picri,ki\u{p2,a2) ^ +M pi,ai) . (5) 

/[p-^ + 1^1 ) 
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Xa in Eq.(l) is an unknown counterterm. We have introduced cutoffs on the invariant mass 
= {p + k)'^ of the two particle sector in the interaction parts of the hamiltonian, Hy and 
Hj^ (see also Q). The integration measure is 

\P) = 4|0) , (7) 
and 6^{p) = S'^{p^)S{p^). Standard light-front commutation relations are 

[ap,al] = 2{27Tfp-^6'{p-k) . (8) 



3 Renormalization 

The similarity transformation Sx transforms -f^A to a band-diagonal hamiltonian Hx, 

Hx = sIHaSx ■ 



(9) 



Expressions for Sx and Hx are found in perturbation theory jT|. X\ in Ha is fitted order by 
order in g, so that Hx does not have A dependent (i.e. divergent) matrix elements for A — s- oo. 
This can be guaranteed in any finite order in perturbation theory. 
In the second order the transformation gives: 



H2 — fx( H^ + X2 - - 



{l-fx)HY,il + fx)HY 



(10) 



The underlining denotes the energy denominator and fx is the diagonal proximum operator (see 
Appendix). 

In the fermion-boson - fermion-boson sector Eq. (10) reads 



H 



-fx[H+- ^ {l-fx)H^ {l + fx)H^ + -(1 + fx)H^ {l-fx)H^ 



[ID 



This expression is not divergent for A ^ cxo, thus no counterterm is needed in this sector. 
However, in the fermion-fermion sector, one obtains 



H 



2/-/ 



-{I- fx)H^H^+X2A . 



(12) 



The loop integration in the first term is linearly divergent. The form of this divergence dictates 
the form of the second order counterterm. Explicitly, one has to choose 



^2A = E/b]M(p^l 



1 9' 



p^ 



167r2 



1 A^ 
-A^ + {?,m'^ - i?)\og^ + A 



(13) 



where A is an undetermined constant. 
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Higher order calculations lead to the following expressions for and X^fh-fbA- 



X3A 




+ E [pi,P2,mA^-Ml^,)2{27rf6'{p,-p2-k)x 



0"l,0'2 




"T- log IT \P2'^2,k)(piai\u{p2,(T2)-—^u{pi,ai) + h.c. 




(14) 





(15) 



where B, C and D are finite unknown constants. 

There is also another term of order in the fermion-fermion part of X^. We did not calculate 
it because our goal was to investigate the possibility of fitting finite parts of counterterms by 
requesting the T-matrix covariance in the fermion-boson channel (Section 4.1) and the emergence 
of the Dirac equation for physical fermions. As A^4/-/a does not contribute either to T4 or the 
second order Dirac equation, it was irrelevant for our considerations. Also, X4j_ja is more 
complicated to calculate than the terms we need to discuss here, because of two correlated loop 
integrations. 

4 Finite parts of the counterterms 

The renormalization procedure does not determine values of the finite parts of counterterms. 
To find them we need to introduce extra conditions. In principle, the constants should be fitted 
to match experiment. It is interesting to look for theoretical requirements of symmetries, which 
may constrain these constants. The T-matrix calculated with the general counterterms ([T^)-([l5[) 
is not automatically covariant. So, the covariance of the T-matrix provides useful conditions. 
Another condition will be provided by requiring that the full Hamiltonian eigenvalue equation 
could be reduced to a free Dirac equation. 

4.1 T-matrix 

We calculate our T-matrix using the formula 



The second order T-matrix has a covariant form and does not depend on the counterterms. 
X\ starts contributing in the fourth order. The explicit A dependence of counterterms cancels 
divergences in the loop integrations in other terms. So, T4 is finite. However, it is not covariant 
automatically. 



T{E) =Hj + Hj 



1 



(16) 



E-Ho + ie 



{p2cr2,k2\Ti\pi(Ti,ki) 




e{A^ - Mi)9{A^ - Ml)2{27rf5^{p2 + k2 - pi - ki)x 
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ri(s)^ + r2(s) + r3(s) 



2(p+ + kt 



u{pi,ai) . 



(17) 



To obtain a covariant result for T4 we demand that the function T^^s) vanishes for arbitrary s. 
Its exphcit form reads 



where 



1 

s — 



1 . C 



m 



{s - m )- log — + 3m log — — A + 16it af{s){s - m ) + 7/(s) 



(18) 



s = {pi + kif = Ml , (19) 

and functions a/(s) and 7/(s) are given in Appendix. As 167r^aj(s)(s — m?) +7/(s) turns out 
to be real and independent of s, the condition T3{s) = implies two relations: 



and 



B = C 



-m + /i log — 5- + 3m lo; 



m^ 



m 



(20) 
(21) 



4.2 Dirac equation 

To describe a physical state in terms of free Fock states one considers the eigenvalue equation 



_|_ 2 

P'^) physical — p_^_ physical 



(22) 



The physical fermion state is a superposition of the bare fermion and fermion-boson states: 

\Pa) physical = J2'^''JP^^) + T. f[p,k]2{27rf5HP-p-k)ra,ix,M^)\pa2,k) . (23) 
By following steps from ref. 0] one can reduce Eq.(p2|) to 



2P^ 



-0 = 



(24) 



for the one-body sector wavefunction ^. Using our hamiltonian with counterterms restricted by 



conditions (|20|)-(21), one gets 



^3 



1 + 



1 + 



9 



167r2 
,2 



167r2 
+ o(/). 



3 A 

-log — -/3(m2) 
a(m^) +o(/) , 



+ o(/), 



(25) 

(26) 
(27) 
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Our earlier demand of the T-matrix covariance established the value of the mass counterterm 
X2 in a way that also leads to the vanishing of H3 in order g^. 

In general, one can expand both non-zero H's in a power series in g 

S = + H(2)<^2 ^ ^{4)^4 ^ . . . (28) 
and one can translate the requirement that m is the mass of physical fermions, 

■m)V = 0, (29) 



^(0) , ^(2) 2 , ^(4) 4 , 

"2 + "2 y I ^2 y ' 



into the condition for all coefficients 

^1 — "2 • 



(30) 

This is the threshold condition which makes the T-matrix threshold to appear at s = {m + /i)^, 
where m is the position of its fermion pole. 

Let us investigate which terms of H contribute to H^*^ If one puts g = then, the only 
condition one gets is 

\P(y)physical = \P(^) ■ (31) 

Technically, the zeroth order terms H^'^^ and Hg'^^ come from the inversion of J2a''^P(^rnUpam = 
"I" wl^ich is a part of H^H^. Dirac equation results in this order automatically; H^*^^ = 
^2 ■ 

"(2) "(2) 

One can easily see that the second order terms and partly come from the term 
/f>_X+3. So, one needs third order vertex corrections, such as Xj^^, to know all second order 
contributions to the Dirac equation. There is an unknown finite parameter D in X4.3. The 
condition = Hg^^ and Eqs. (p5|)-(p6[) lead to 

log^ = ^-167r2[a/(m2)+/3j(m2)] . (32) 

The functions a/(s) and /3/(s) are given in Appendix. 

We see that the requirement that m is equal to the mass of physical fermions implies one 
more condition on the free parts of counterterms. 



4.3 Discussion 



Collecting conditions (^), ( |2l| ) and (32) together, and looking at the structure of the countert- 
erms, we can observe the following. can be accounted for by changing the coupling constant 



9^9 + 
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A2 



in the original hamiltonian, while -'^4/fe-/feA shifts g^ in the seagull term Hj^: 

.4 ^2 



9 



9'+ ' 



32.2 log^^ 



(33) 



(34) 



6 



So, these two counterterms can be absorbed in one, A-dependent coupling constant (|3^) . We 
have to stress that, in physical results, A dependent logarithms log cancel out, leaving 



Thus, g and C will never appear independently, and we have one parameter, combination (35), 
that can be fixed from experiment. 

X2 shifts the mass in the one fermion free energy. Sum of Hy and one of the third order 
counterterms, X^^, reproduces the same uu coupling but with shifted mass of the spinor in the 
one-particle sector, according to the formula 

I 1 + j Um{p, Cr) = Um+Sm{P, O") • (36) 



5 Comparison with Ref. [3] 

It was shown in Ref. Q that, in this model, to get finite and covariant results for the T-matrix 
to all orders of perturbation theory, and to get the mass in the Dirac equation which is required 
by the threshold condition, it is enough to (1) add to the bare cut-off hamiltonian a term that 
shifts the mass of fermions in the free part Hqj, (2) correspondingly, change the spinor mass in 
the vertex, see Eq. (^6|), and (3) allow the coupling to depend on A. 

When one rewrites the hamiltonian of Ref. Q using the invariant mass cutoff and expands it 
in powers of g{m?) up to the fourth order, one gets the same result as obtained in our similarity 
calculation with 

replaced by 

5("i) - -jf{m)af{m^) . (38) 
So, one can choose C leading to the same result as in Ref. [^. 



6 Conclusion 

This work provides an example of application of the similarity renormalization scheme in its 
algebraical version. We have shown how this systematic procedure leads from a divergent hamil- 
tonian to a finite one. The finite hamiltonian gives a covariant scattering matrix in perturbation 
theory. 

The hamiltonian we used was known to lead to covariant results when one introduced special 
counterterms. The question was if a systematic procedure, the similarity renormalization scheme, 
would produce the same solution. The answer is yes. 

On the other hand, Ref. [Q] has recently suggested that the model may find applications in 
pion-nucleon physics when another Fock sector, with one fermion and two bosons, is included. 
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Therefore, our work also suggests that a systematic improvement in the hght-front hamiltonian 
approach to relativistic nuclear physics may be achivable using the similarity renormalization 
group techniques. 
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Appendix 

For any operator A, 

A = l |l)(2|^i2, (39) 

A is defined as 

A= [ |1)(2|— J— ^12, (40) 
where £"s are eigenvalue of Hq. A is a solution of an equation: 

[A,Ho]=A. (41) 
Action of diagonal proximum operator fx is defined as follows : 

hA = j |1)(2|/a(1,2)Ai2. (42) 

We have chosen 

h{l,2)=e{\'-\Ml-Ml\). (43) 



Functions a(s), I3{s) and 7(5) are defined by 

a{s) = ^ / dM^ dx e{K^ - M'^ ^ 

lOTT^ J 



pfs) = --1- / dM^ dx eiK^ - M^) , 

'^^ ' 167r2 j ^ > -s + ie' 

2 ,,2^ {l-x)M^-i? + {l-x)m'' 



7(s) = j dM^ dx 9{k^ - 



- s + ie 



where x is integrated over the whole kinematically allowed region. Their finite parts are defined 

by 



a/(s) = lim 

A— >CXD 
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{s — m^- 2ii^) log —2 



m- 
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